Observational constraints on assisted k-inflation 
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We study observational constraints on the assisted k-inflation models in which multiple scalar 
fields join an attractor characterized by an effective single field (f>. This effective single-field system 
is described by the Lagrangian P = Xg(Y), where X is the kinetic energy of (j>, A is a constant, 
and g is an arbitrary function in terms of Y = Xe x ^. Our analysis covers a wide variety of 
k-inflation models such as dilatonic ghost condensate, Dirac-Born-Infeld (DBI) field, tachyon, as 
well as the canonical field with an exponential potential. We place observational bounds on the 
parameters of each model from the WMAP 7yr data combined with Baryon Acoustic Oscillations 
(BAO) and the Hubble constant measurement. Using the observational constraints of the equilateral 
non-Gaussianity parameter /^l" > we further restrict the allowed parameter space of dilatonic ghost 
condensate and DBI models. We extend the analysis to more general models with several different 
choices of g(Y) and show that the models such as g(Y) = Co + c p Y p (p > 3) are excluded by the 
joint data analysis of the scalar/tensor spectra and primordial non-Gaussianities. 



I. INTRODUCTION 



The cosmic acceleration in the early Universe- 
inflation- has been the backbone of the high-energy cos- 
mology over the past 30 years. In addition to addressing 
the horizon and flatness problems plagued in Big Bang 
cosmology inflation generally predicts almost scale- 
invariant adiabatic density perturbations Q (see for 
reviews). This prediction is consistent with the observa- 
tions of the Cosmic Microwave Background (CMB) tem- 
erature anisotropies measured by COBE Q and WMAP 
It is possible to distinguish between a host of infla- 
tionary models by comparing the theoretical prediction 
of the spectral index n s of curvature perturbations and 
the tensor-to-scalar ratio r with observations, but still 
the current observations are not sufficient to identify the 
best model of inflation. 

In the next few years, the measurement of CMB tem- 
perature anisotropies by the PLANCK satellite Q will 
bring more high-precision data. In addition to the possi- 
ble reduction of the tensor-to-scalar ratio to the order of 
0.01, the non-linear parameter /nl of primordial scalar 
non-Gaussianities may be constrained by about one or- 
der of magnitude better than the bounds constrained by 
the WMAP group. This can potentially provide further 
important information to discriminate between many in- 
flation models. 

The conventional single-field inflation driven by a 
canonical scalar field with a potential V(4>) predicts 
small primordial non-Gaussianities with /nl of the order 
of slow-roll parameters (see [l2[ for early works) . 

However the kinetically driven inflation models (dubbed 
"k-inflation" [l3|) described by the Lagrangian density 
P(4>,X), where X is the field kinetic energy, can giv e 
rise to large non-Gaussianities with |/nl| 3* 1 [lil [l5j. 
This is related to the fact that for the Lagrangian in- 
cluding a non-linear kinetic term of X the propagation 
speed c s is different from 1 (in the unit where the speed 



of light c is 1) [16H18J- Since the non- linear parameter is 
approximately given by /nl ~ — 1 jcr s > one nas I /nl | 3> 1 
for c 2 < 1. 

In the models motivated by particle physics such as 
superstring and supcrgravity theories, there are many 
scalar fields that can be responsible for inflation [H, Q . 
In some cases, even if each field is unable to lead to 
cosmic acceleration, the presence of many fields allows 
a possibility for the realization of inflation through the 
so-called assisted inflation mechanism fl9j . In fact, mul- 
tiple (canonical) scalar fields with exponential potentials 
Vi((f>i) = Cje ^ evolve to give dynamics matching a 

sing le field with the effective slope A = (J2i=i V^f) 
[19| . Since A is smaller than the individual Aj, the pres- 
ence of multiple fields can lead to sufficient amount of 
inflation H^. 

If we take into account a barotropic perfect fluid (den- 
sity p m ) in addition to the canonical scalar field (den- 
sity with the exponential potential V{(f>) = ce~ A *, 
there exists a so-called scaling solution along which the 
ratio p^/pm is constant [2ll. |22|. In the presence of non- 
rclativistic matter the scaling solution is unstable for 
A 2 < 3, in which case another scalar-field dominated 
solution is a stable attractor (22[. If A 2 < 2, the lat- 
ter can be used for inflation as well as dark energy. If 
we extend the analysis to the models described by the 
general Lagrangian P{4>, X) then the condition for the 
existence of scaling solutions restricts the form of the 
Lagrangian to be P = Xg(Y), where A is a constant 
and g is an arbitrary function in terms of Y = Xe x ^ 
[H Provided A 2 < 2dP/dX there exists a scalar- 
field dominated attractor that can be responsible for in- 
flation pol [26j . In fact this Lagrangian covers a wide 
class of inflationary models such as the canonical scalar 
field with the exponential potential (g(Y) = 1 — c/Y, i.e. 
P = X — ce~ x ^) [27j and the dilatonic ghost condensate 
model MY) = -1+cY, i.e. P = -X + ce^X 2 ) ||| (see 
Refs. [28|, [29| for the original ghost condensate model). 
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In the presence of multiple scalar fields it was shown 
that the Lagrangian P = Y17=i -^idO^i)) where giYi) is 
an arbitrary function with respect to Yi = JQe ^ , gives 
rise to assisted inflation [25|, |3(| , as it happens for the 
canonical field with the exponential potential. In other 
words, in the regime where the solutions approach the as- 
sisted inflationary attractor, the system can be described 
by the effective single-field Lagrangian P = Xg(Y) with 

Y = Xe^ and the slope A = (£" =1 1/A 2 )~ 1/2 . While 
the scalar propagation speed is different from I in those 
models, the scalar spectral index n s and the tensor-to- 
scalar r are written in terms of the function g(Y) and 
its derivatives g'(Y), g"(Y). By specifying the func- 
tional form of g(Y), the observablcs n s and r as well 
as the equilateral non-Gaussianity parameter /^l" can 
be expressed by the single parameter A in the attractor 
regime. This property is useful to place tight observa- 
tional bounds on those models. 

In this paper we confront the assisted k-inflation sce- 
nario described by the effective single-field Lagrangian 
P = Xg(Xe x ^) with the recent CMB observations by 
WMAP 0] combined with BAO HJ and the Hubble 
constant measurement (HST) [32| • We evaluate three ob- 
servables n a , r, and /nl" without specifying the forms of 
g(Y) and apply those results to concrete models of infla- 
tion. We place observational constraints on a number of 
assisted inflation models such as (A) canonical field with 
the exponential potential, (B) tachv on 1 3311 . (C) dilatonic 
ghost condensate, and (D) DBI field [34[. Since the effect 
of the non-linear term in X is important in the models 
(C) and (D), the primordial non-Gaussianity can reduce 
the parameter space constrained by the information of n s 
and r. 

We shall also study other assisted inflation models such 
as g(Y) = co + J2 P ^o c v^ v an d the generalization of the 
DBI model. Interestingly the observational bound from 
the equilateral non-Gaussianity parameter f^* 1 com- 
bined with n s and r can rule out some of those models. 



II. BACKGROUND DYNAMICS IN ASSISTED 
K-INFLATION 

We start with the single-field k-inflation models de- 
scribed by the action fl3l j 



S 



(Pxy/-g M 



R 
~2 



PU,X) 



(1) 



where gu is a determinant of the metric c/ M „, R is a scalar 
curvature, P is a general function in terms of the scalar 
field 4> and the kinetic term X = —g^d ll (f>d^<f>/2. We use 
the unit M pl = 1, where M p \ = (8-kG)" 1 / 2 is the reduced 
Planck mass (G is gravitational constant), but we restore 
Mpi when the discussion becomes more transparent. 

The pressure P and the energy density p of the field (f> 
are given, respectively, by 



where P.x = dP/dX. We also define the equation 
of state 'wj,, as = P/p = P/{2XP iX - P). The 
cosmic acceleration can be realized under the condition 
\2XP X \ < \P\, i-e. either (i) X is small, or (ii) P x is 
small. The case (i) corresponds to conventional slow-roll 
inflation driven by a field potential, whereas the case (ii) 
to kinetically driven inflation fl3j ] . One of the examples 
in the class (ii) is the ghost condensate model [111 [29[ de- 
scribed by the Lagrangian P = —X + X 2 /M A , in which 
case inflation occurs around X = M i /2. 

In Refs. [23|,[24| it was shown that the condition for the 
existence of cosmological scaling solutions in the presence 
of non-relativistic matter restricts the Lagrangian of the 
form 



P(^,X)=Xg(Y), Y = Xe 



(3) 



where A is a constant and g is an arbitrary function in 
terms of Y. This Lagrangian was derived by imposing 
that Q<p I f2 m = constant and = constant in the scaling 
regime (where VL^ and VL m are the density parameters of 
the scalar field and non-relativistic matter, respectively). 

For the Lagrangian © there is another solution that 
can be responsible for the cosmic acceleration. This cor- 
responds to the fixed point with the equation of state 

m 



-1 



A 



3P 



x 



(4) 



The condition for the cosmic acceleration is w,p < —1/3, 
i.e. A 2 < 2P t x- Since this point is stable for A 2 < 3Px, 
the solutions approach it provided that inflation occurs. 
Under the condition A 2 < 3P.x the scaling solution is 
unstable [25j . 

Let us consider the models with multiple scalar fields 
4>i (i = 1, 2, • • • , n) described by the action 



S = Id xy/-g M 



E 



(5) 



where Xj = —g^d tl (j>id u <j>i/2, A^'s are constants, and g is 
an arbitrary function in terms of Yi = Xie x '^' i . Since we 
focus on inflation in the early Universe, we do not take 
into account other matter sources in the action ([5]). In 
the flat Fricdmann-Lcmaitrc-Robertson- Walker (FLRW) 
background with a scale factor a(i), the equations of mo- 
tion are 



3H 2 



i=l 



2H=-J2{Pi + Pi), 

i=l 

p i + 3H(p i + P i ) = 0, 



(6) 



(7) 



(i = l,2,...,n), (8) 



where H = a /a is the Hubble parameter (a dot denotes 
a derivative with respect to t), and 



P = P(<j>,X), 



2XP 



x 



P, 



(2) 



Pi = Xig(Yi), Pi = X, [g{Yi) + 2Y i g'(Y i )] 



(9) 
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Here and in the following a prime represents a derivative 
of the corresponding quantities, e.g., g'{Yi) = dg/dYi. 

In order to discuss the cosmological dynamics for the 
theories described by the action ([5]) we introduce the fol- 
lowing quantities 



then Eq. (|2TJ yields 



V3H 



(10) 



The differential equations for the variables Xi and j/j are 
given by 



dxi 

In 



~2 



" a/6 



V6A(Yi) 



dyi % 
dN ~ 2 



V6, 



(11) 
(12) 



where A" = In a is the number of c-foldings, and 

MY) = [g(Y i ) + W i g'(Y i )+2Y?g"(Y i )]-\ (13) 
O 0j = z? [ 9 (Y) + 2Y i5 '(Y i )] . (14) 

From Eqs. ([TT]) and fj 12[) we find that the fixed point 
{dXi/dN = and dyi/dN = 0) responsible for inflation 
(z/j 7^ 0) satisfies 



A; 



g(Y i )+2Y i g'(Y i ) 2 



1 + ^(7(^)3 



Then the equation of state for each field, w l p i 
g(Y i )/[g(Y i ) + 2Y l g'(Y i )}, reads 



(15) 



w <t>i 



T + — MXi 



(16) 



We require that Eq. (fT5"j) is satisfied for all i = 1, 2, • • ■ , n. 
Hence A independent of i, i.e. 



X 1 x 1 = ■ ■ ■ = \ l x l = • • • = \ n x n = Xx . 
This property also holds for Yi and w c j >i : 
Yi = ■ ■ ■ = Yi = ■ ■ ■ = Y n = Y , 

Wfa — ■ ■ ■ = Wfc = ■■■ = W<f> n = . 

From Eq. (fl"5j) it follows that 

V6[g(Y) + Yg'(Y)] 



Xx = 



g(Y) + 2Yg'(Y) 
V6 



If we choose 



1 

l+.g(y)x 2 A 2 ^ ' 



1 - 1 

A 2 = ^ A? ' 



(17) 

(18) 
(19) 

(20) 
(21) 

(22) 



Xx = ^[l + g(Y)x 2 } 



(23) 



This shows that, along the inflationary fixed point, 
the system effectively reduces to that of the single 
field with the Lagrangian P = Xg(Y) with Y — 
Xe x< ^ . Since the sum of the density parameters = 
x 2 [g(Yi) + 2Yig'(Yi)] satisfies the relation = 1' 

we have 



x 2 [ g (Y) + 2Y g '(Y)} = i 



(24) 



From Eq. (gU]) it then follows that x = X/(y/6P x ), where 
we have used P t x = g(Y) + Yg'(Y). The field equation 
of state (fT6|) is given by 



W0 = -1 + 



A 



3P 



x 



(25) 



From Eq. (f2"2"j) we find that the effective slope squared 
A 2 is smaller than A 2 of each field. Even when the cos- 
mic acceleration does not occur with a single field, it is 
possible to realize inflation in the presence of multiple 
fields. The above discussion shows that assisted inflation 
occurs for the multi-field k-inflation models described by 
the action (JSJ) - In the regime where the solutions ap- 
proach the assisted inflationary attractor satisfying the 
condition A 2 < 2Px, the multi- field system reduces to 
that of the effective single field. In the following we shall 
study the effective single-field system described by the 
Lagrangian with the slope A given in Eq. (|2"2"j) . As 
we have mentioned in Introduction, this analysis covers 
a wide variety of assisted inflation models. 



III. INFLATIONARY OBSERVABLES 

It is possible to distinguish between a host of infla- 
tionary models by considering the spectra of primordial 
density perturbations generated during inflation. For the 
calculations including primordial non-Gaussianitics it is 
convenient to use the ADM metric (3|| of the form 

ds 2 = - [(1 + a) 2 - a- 2 {t)e~ 2n {d%lj) 2 ] dt 2 + 2d i ipdtdx i 
+a 2 (t)(e 27 %- + h ij )dx i dx j , (26) 

where a, ip, and 1Z are scalar perturbations, and fay are 
tensor perturbations. We do not take into account vec- 
tor perturbations because they rapidly decay during in- 
flation. 

In the metric ([26)1 we have gauged away a field E that 
appears as a form E ij inside the last parenthesis. This 
fixes the spatial part of the gauge-transformation vector 
We also choose the uniform- field gauge such that the 
inflaton fluctuation Scf> vanishes (5<f> = 0), which fixes the 
time component of £ M . 

Integrating the action ([T]) by parts for the metric 
(|26p and using the background equations of motion, the 
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second-order action for the curvature perturbation can 
be written as nil 



So 



dt d 3 x a 3 Q 



where Q = e/cr s , and 
H 



P 



x 



P 



X 



2XP 



XX 



(27) 



(28) 



The conditions for the avoidance of ghosts and Laplacian 
instabilities correspond to Q > and c 2 > 0, respectively, 
which are equivalent to 



e > and 



ct > 0. 



(29) 



For the Lagrangian including a non-linear term in X (i.e. 
P.xx 7^ 0), the scalar propagation speed c s is different 
from 1. 

The equation for the Fourier mode of 1Z follows from 
the action (f2"7|) . For the modes deep inside the Hubble 
radius we choose the integration constants of the solution 
of 1Z to recover the Bunch-Davies vacuum state. After 
the perturbations leave the Hubble radius (c s k < aH, 
where A; is a wave number) the curvature perturbation is 
frozen, so that the scalar power spectrum is given by 



_L IP 

8^ 2 M 2 , ^7 



(30) 



which is evaluated at c s k = aH. The scalar spectral 
index is 



n s - 1 = 



dlnK 



where 



din A: 



e 
lie 



c s k—aH 



-2e-r)-s, (31) 



(32) 



He, 



Here we have assumed that the field propagation speed 
slowly changes in time, such that \s\ <^ 1. 

For the theories described by the action ([1]) the tensor 
perturbation hij satisfies the same equation of motion as 
that for a massless scalar field. Taking into account two 
polarization states, the spectrum of hij and its spectral 
index are given, respectively, by fl6j 



2H 2 
dlnVt 



dink 

The tensor-to-scalar ratio is 



= -2e. 



k=aH 



16c s e 



ic s n t . 



(33) 
(34) 

(35) 



The non-Gaussianity of the curvature perturbation is 
known by evaluating the vacuum expectation value of 
the three-point correlation function (TZ(ki)TZ(k2)TZ(k3)) , 
where 7Z(ki) is the Fourier mode with a wave num- 
ber kj (i = 1,2,3). We write the bispectrum 
in the form (n(k 1 )n{k 2 )TZ(k 3 )} = (2Tr) 3 S^(k 1 + 
k 2 + k 3 ){V s ) 2 B(k 1 ,k 2 ,k 3 ), where k t = \k t \. In k- 
inflation one can take a factorizable shape function 
B = (2^) 4 (9/ NL /10)hl/(fc?fc 2 3 ) - l/(k 3 k 3 ) - l/(k 3 k 3 ) - 
2/{k\k 2 k 2 ) + l/(fcifc|fc|) + (5 perm.)], where the permu- 
tations act on the last term in parenthesis (36l . l37j | . For 
the equilateral triangl es ( k\ = k 2 = k 3 ), the non-linear 
parameter is given by [lj, [HI, [l8j 



/■cquil 



5_ f }__ _2fJ. 

1 ,2 g 



81 V c- 
55 e 
+ 36cl 



5 7] 

12g 



35 
" 108 

85 s 
54 c|' 



1 

-77-1 



where 



E 
/' 



XP X + 2X 2 P, XX = H 2 e/c 2 s , 

X 2 Pxx + 2X 3 P, XX x/3 
S / 1 2 e s 



G 



1 



3 e x c 2 



(36) 

(37) 
(38) 



and e x = -(X/H 2 )(dH/dX). In the second line of 
Eq. (|38| we have used X = —QHc 2 Xex/e, which follows 



from the background equation of the field (f> |14| . Our sign 



convention of /nl" 1 ' coincides with that in the WMAP 
7yr paper 0. The observational bound on the equilat- 
eral non-linear parameter constrained by the WMAP 7yr 
data is 



214 < /° q L uil < 266 (95%CL). 



(39) 



Let us consider the case in which the multiple fields 
join the effective single-field attractor characterized by 
the conditions (|2"0 ]) -([2l[ ]) . From Eqs. (|23J) and (JHJ) we 
obtain 



A 



6[g(Y) + Yg>(Y)}< 
g(Y) + 2Yg'(Y) 



(40) 



By choosing a specific function g(Y) and solving Eq. (|4"0"j) . 
we can determine Y in terms of A (i.e. Y is constant). 
The slow-roll parameter e and the scalar propagation 
speed squared c 2 are 



3[g(Y) + Yg'(Y)] 



g(Y) 



■2Yg'(Y) ' 
j(Y) + Yg'(Y) 



g(Y)+5Yg'(Y) + 2Y 2 g"(Y) 



(41) 
(42) 



which arc functions of Y only. Then one has e = constant 
and c 2 = constant on the inflationary attractor, thereby 
leading to r\ = 0, s = 0, and /x/S = (1/c 2 — l)/6. From 
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Eqs. (J3TD, (041), ([35 

servables reduce to 

n s - 1 = 

r = 



and ([36)1 the three inflationary ob- 



-2e = 
16c s e 

/•cquil 275 



n t . 



8c s (l-n B ) ; 



1-1 



55 e 
36 cl 



(43) 
(44) 

(45) 



where e and c 2 are given in Eqs. (|4Tj) and (|42]) . Since e 
is constant along the inflationary attractor, there are no 
runnings for scalar and tensor perturbations. 

Since Y is known in terms of A for given g(Y), all 
the observables in Eqs. (|43l) - (|45"|) can be expressed by A 
(or e). Observationally one can place the bounds on the 
parameter A for each model. In the following we shall 
proceed to the observational constraints on assisted k- 
inflation models. 



IV. OBSERVATIONAL CONSTRAINTS ON 
FOUR MODELS OF ASSISTED INFLATION 

In this section we study the observational constraints 
on a number of assisted inflation models by choosing spe- 
cific forms of g(Y). These include (A) canonical field 
with an exponential potential [g(Y) = 1 — c/Y], (B) 
tachyon [g(Y) = -cy/1 - 2Y/Y], (C) dilatonic ghost 
con densate [g(Y) = -1 + cY], and (D) DBI field [g(Y) = 
— — 2Y /Y — c/Y], where c is constant. 

In the model (A) one has (? s = 1, so that the non- 
Gaussianity is small enough (f^ d — 55e/36 <C 1) to 
satisfy the observational bound (|39|) . We can constrain 
either A or e by carrying out the CMB likelihood analysis 
with respect to n s , r, and n t . In the tachyon model (B) 
the scalar propagation speed c s does not equal to 1, but 
the difference from 1 is required to be small. Hence the 
situation is similar to that in the model (A). 

For the models (C) and (D) c s can be much smaller 
than 1 , while satisfying the condition e«l. In such cases 
it is possible to place tight bounds on the models from the 
primordial non-Gaussianitics in addition to those coming 
from n s , r, and nt. 



A. Canonical field with an exponential potential 



The canonical field with the exponential potential de- 



scribed by the Lagrangian P 
to the choice 



X 



,-\<t> 



g(Y) = l-c/Y. 



corresponds 



(46) 



In this case one has c/Y = 6/A 2 — 1, e = A 2 /2, and 
c 2 = 1. Inflation occurs for A 2 <C 1, i.e. X < ce~ A *. 



The inflationary observables are 

n s — 1 = n t = — A 2 

r = 8A 2 , 
fZ n = 55A 2 /72. 



(47) 
(48) 
(49) 




Figure 1: One-dimensional marginalized probability distribu- 
tion of the parameter A for the model P = X — ce~ A ^ con- 
strained by the joint data analysis of WMAP 7yr, BAO, and 
HST. We use the theoretical expression of n s , r, and n t given 
in Eqs. |07| and 



Using the Cosmological Monte Carlo (CosmoMC) code 
[38[, we carry out the likelihood analysis with the 
WMAP 7yr data combined with BAO and HST. As we 
show in Fig. [TJ the likelihood analysis in terms of n s , nt, 
and r gives the following bound 



0.086 < A < 0.228 



(95% CL). 



(50) 



The Harrison-Zcl'dovich (HZ) spectrum (n s = 1 and r = 
0) is disfavored from the data. Under the bound (f5Tj|) one 
has /nl" 1 ' ^ 1) sucn that the non-Gaussianity constraint 
is satisfied. 



B. Tachyon 

A tachyon field ip with a potential V((p) corresponds 
to the Lagrangian P = — V(f)yl 



-g^d^d^/2 



2X, where X 
Choosing the function 



g(Y) = -cVT^W/Y, 



(51) 



where Y = Ae A * = X, one can show that the 
Lagrangian P = Xg(Y) reduces to the form P = 

— 4c/(A 2 </3 2 )\/ 1 — 2X. Hence the tachyon potential 
V(ip) oc ip~ 2 leads to assisted inflation. The cosmolog- 
ical dynamics in the presence of the inverse power-law 
tachyon potential have been discussed in Refs. |39j . 

= 3Y a nd c 2 = 
6cY/\/l-2Y. 



For the choice (|5T|) it follows that e 
1 — 2Y, where Y is related to A via A 2 = 
The inflationary observables are 



-6F, 



(52) 
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r = ASYy/l 

f cquil _ 3905 _ 



2Y. 



Y 



972 1 -2Y 



(53) 
(54) 



Since we require 7 C 1 to realize the nearly scale- 
invariant scalar spectrum, the non-Gaussianity is sup- 

<C 1). The relation between 
1 - n„Vl - 2Y ~ 8(1 



cquil 



pressed to be small (/ NL 
r and n s is given by r = 
which, in the limit that Y 



n s )vT-2F ~ 8^1 - n, 
0, is the same as that for the 
canonical field with the exponential potential. This prop- 
erty comes from the fact that tachyon inflation is driven 
by the potential energy rather than the field kinetic en- 
ergy. The joint CMB likelihood analysis combined with 
BAO and HST gives the bound 

1.7 x 1CT 3 < Y < 7.7 x 1CT 3 (95% CL). (55) 

Then c 2 = 1 — 2Y is indeed close to 1. 




C. Dilatonic ghost condensate 

The dilatonic ghost condensate model is described by 
the Lagrangian P = —X + ce^X 2 , i.e. 



g(Y) = -l + cY. 

In this case we have 

_ 3(2cY - 1) 2 _ 2cY - 1 

6 ~~ 3cY - 1 ' ° s ~ 6cY- 1 ' 

where cY is known by solving Eq. (|40j) . i.e. 



2cF-l = /(A), /(A) = - 



A4 + y A2 



(56) 



(57) 



(58) 



In Eq. ([55]) we have chosen the solution with cY > 1/2 
to avoid the appearance of ghosts (23|. The inflationary 
observables are given by 



1 = nt 



A 2 



/(A)' 



8A 2 



v //(A)[3/(A) + 2]' 



rcquil 



275 
486 



/(A) 



55 A 2 
72 /(A) 



/(A) 



(59) 



(60) 



(61) 



In the limit that A -> one has /(A) ~ A/2^ 
and hence /nl"' — * — 00 ■ Using the WMAP 7yr bound 
/nl UU > -214, we obtain the constraint A > 8.4 x 10~ 3 
(95 % CL). 

In the region A 2 <C 1 one has /(A) ~ V3X/6, n s ~ 
1 - 2^3A, and r ~ 8 • 3 1/4 A 3/2 , which give the relation 
r ~ (2 v / 6/3)(l - n s ) 3 / 2 . In this model the tensor-to- 
scalar ratio is smaller than the order of 0.1, so that the 
allowed region of A is mainly determined by n s . The 
CMB likelihood analysis in terms of n s , n t , r shows that 



Figure 2: One-dimensional marginalized probability distri- 
bution of the parameter A in the dilatonic ghost condensate 
model constrained by the joint data analysis of WMAP 7yr, 
BAO, and HST. We also show the bound on A coming from 
the WMAP 7yr constraint of the equilateral non-Gaussianity 
parameter, S > —214, as well as the bound on A corre- 



sponding to the constraint f\ 



cquil 



> -100. 



A is constrained to be 4.0 x 10~ 3 < A < 1.5 x 10~ 2 
(95 % CL), sec Fig. [5J Combining this with the non- 
Gaussianity constraint, it follows that 



.4 x 10~ 3 < A < 1.5 x 10" 



(95% CL). (62) 



If the future observations constrain the non- 



Gaussianity parameter at the level /, 



cquil 
NL 



> 



-100, it 



will be possible to exclude the dilatonic ghost condensate 
model (see Fig. Moreover the precise measurement of 
the scalar index n s can reduce the allowed range of A 
further. 



D. DBI field 

The DBI field <fi is characterized by the Lagrangian 

P = -/(^)- Vi - W)X + /WO" 1 - V{4>) , (63) 
where f(<p) and V(4>) are functions of (f>. If we choose 

g {Y) = -VI - 2Y/Y - c/Y , (64) 

the Lagrangian P = Xg(Y) reduces to (|63|) with f((f>) = 
e x< t> and V(<f>) = (c+ l)e- A< ^. Hence the DBI field with 
the exponential potential V(4>) = (c + l)e~ A< ^ leads to 
assisted inflation. 

For the function (f6"4l it follows that 



3Y 



cVl -2Y+1 



1-2Y . 



(65) 
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If c < 1 , one has £«1 and c 2 ~ 1 for Y <C 1 • This case is 
similar to tachyon inflation in which cosmic acceleration 
is driven by the field potential. One can also realize e«l 
under the following condition 



VI - 2Y > 1 



(66) 



If c ^ 1, then it is possible to satisfy (|66|) even for 
the values of Y close to 1/2. In fact this is the ultra- 
relativistic regime of the DBI inflation in which the 7 

factor 7 = 1/ \Jl — f{<p)4> 2 is much larger than 1 . Even 
in this "fast-roll" regime the presence of the potential is 
important to satisfy the condition (|66p . 

The inflationary obscrvablcs are 



"1 



/•cquil 
JNL 



3(1 -eg) 
cc s + l ' 
24c s (l - a 
cc s + l ' 
55 (10cc s -71) 
~1944 (cc s + l) 



- 1 



(67) 
(68) 
(69) 



where c s = y/l — 2Y. These obscrvablcs depend not only 
on c s (or Y) but on the coefficient c associated with the 
field potential. For larger c it is possible to satisfy the 
observational constraints of n s , r, and n t with smaller 
c s , because the denominators of Eqs. ([6T)) and ([6"8")l get 
larger. In fact, Fig. [3] shows that, for larger c, the one- 
dimensional marginalized probability distribution of A 
tends to shift to the regions of smaller c s . In Fig. [3] the 
propagation speed c s close to 1 is not favored because n s 
and r are close to the HZ spectrum. The models with 
very small c s are also disfavored because of the large de- 
viation from the HZ spectrum. 

In Fig. U we plot the non-Gaussianity parameter fS^ 1 
given in Eq. (|69|) versus the scalar propagation speed c s 
for three different values of c. For c = 10 2 we obtain the 
bound c s > 3.2 x 10~ 2 from the WMAP 7yr upper limit 
/n1 u " < 266. On the other hand, the WMAP 7yr lower 
limit 11 > -214 gives the bounds c s > 3.1 x 1CT 2 and 
c s > 3.6 x 10~ 2 for c = 10 3 and c = 10 4 , respectively. 

As we see in Fig. [3l the CMB likelihood analysis in 
terms of n s , r, and n t places the constraints on c s , 
as 0.48 < c s < 0.84 (95 % CL) for c = 10 2 and 
0.06 <c s < 0.35 (95 % CL) for c = 10 3 . If c < 10 3 the 
non-Gaussianity does not provide additional constraints 
on c s to those derived by the likelihood analysis in Fig. [3] 
If c > 10 3 the non-Gaussianity plays an important role to 
restrict the allowed parameter space of c s further. In par- 
ticular, for c = 10 4 , there are almost no allowed regions 
to satisfy all the observational constraints (see Fig. [3]). 
Hence the models with c > 10 4 are excluded by the anal- 
ysis including non-Gaussianitics. 




Figure 3: One-dimensional marginalized probability distribu- 
tion of the field propagation speed c s (with the logarithmic 
scale) in the DBI model constrained by the observational data 
of WMAP 7yr, BAO, and HST. The three solid lines corre- 
spond to the cases: (i) c = 10 2 , (ii) c = 10 3 , and (iii) c = 10 4 . 
We also show the bound derived from the non-Gaussianity 
constraint > —214 in the limit cc s 2> 1. 








0.05 



Figure 4: The equilateral non-Gaussianity parameter f^ 11 
versus the scalar propagation speed c 3 in the DBI model for 
(i) c = 10 2 , (ii) c = 10 3 , and (iii) c = 10 4 . For c = 10 2 the 
scalar propagation speed is constrained by the WMAP 7yr 
upper bound < 266, whereas for c = 10 3 and c = 10 4 it 



is constrained by the lower bound f ] 



cquil 
NL 



> -214. 
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V. MORE GENERAL MODELS 

So far we have studied the observational constraints on 
four assisted inflation models. Among them the dilatonic 
ghost condensate and the DBI models can be tightly con- 
strained by taking into account the bound coming from 
the primordial non-Gaussianity. This is associated with 
the fact that both e and c 2 can be much smaller than 1 in 
those models. In this section we shall extend the analysis 
to more general functions of g(Y) . 

In the dilatonic ghost condensate model the numera- 
tors of e and c 2 in Eq. (|57[) vanish at cY = 1/2, whereas 
the denominators of them are non-zero finite values. In 
the DBI model the numerator of e in Eq. ((65)) does not 
vanish in the ultra-relativistic regime (Y rj 1/2), whereas 
c 2 <C 1. In the DBI case it is possible to have 6 C 1 as 
long as the denominator of e is much larger than the 
numerator of it [which is satisfied under the condition 
(|6"6"]) ]. Since these models are qualitatively different, we 
classify the assisted k-inflation models into two classes in 
the following discussion. 



A. Class (i) 

Let us first study the models in which inflation occurs 
around Y = Yq, where Yq satisfies 



g(Y o ) + Y o g'(Y o ) = 0. 



(70) 



As in the case of the dilatonic ghost condensate, we con- 
sider the models in which the numerators of e and c 2 
in Eqs. (|4Tj) and ([42]) vanish, whereas the denominators 
are non-zero. Since Y = Yq corresponds to the exact de 
Sitter solution, we perform the linear expansion of the 
variables e(Y) and c 2 {Y) by setting Y ~ Yq + SY with 
\SY/Y \ < 1. It then follows that 



e(Y) 



e'(Y )SY=^ 



1 



Y g"(Y ) 



WW) 



c 2 s (Y) ~ c 2 s '(Y )SY 



2Yn 



■SY. 



SY, (71) 
(72) 



This shows that the ratio c 2 /e is approximately constant 
in the regime |5F/F | < 1: 



1 

12 



Y g"(Y ) 



Expanding Eq. ([40]) at Y = Yq, we have 
Y g'(Y ) 



(8YY = 



6{2g>(Y )+Y g"(Y )}< 



A 



(73) 



(74) 



As long as g'(Yo) > there exists a solution with SY > 
0. The conditions (f2T)]) for the avoidance of ghosts and 
Laplacian instabilities translate into 



g'(Y ) > o, 

Y g"(Y ) > -2g'(Y ) 



(75) 
(76) 



In the ghost condensate model described by the func- 
tion giY) = — 1 + cY the second derivative g"(Y) au- 
tomatically vanishes, which gives cf/e = 1/12. In this 
model the variable A is observationally bounded as A < 
1.5 x 10~ 2 (95 % CL), in which case SY/Y Q = A/(2 V / 3) < 
4.3 x 10 -3 . Hence it is a good approximation to use the 
linear expansion given above. In fact we have carried out 
the CMB likelihood analysis by employing the relation 
c 2 /e = 1/12 and confirmed that the observational bound 
on A is very similar to that given in Eq. (1621) . 

We study the following more general models 



g(Y)=c + J2> 

p#0 



(77) 



where c p are constants. The power p can be integer or 
some real number. From Eq. (|73[) the ratio c 2 /e is given 
by 



=2 i Ep^r 1 



(78) 



For the single power p, i.e. giY) = c + c p Y p , Eq. ([78 
reduces to 



1 



6(p + l) 



(79) 



The conditions (|75|) and f|T6[) give pc p > and p(p + 
l)c p > 0, respectively, which demand that p > — 1. More 
precisely we require c p < for — 1 < p < and c p > 
for p > 0. In Fig. [5] we plot the line ([79]) in the (e,c 2 ) 
plane for five different values ofp (= —0.5,0.5,1,2,3). 
The ghost condensate model corresponds to p = 1 with 
the tangent c?Je = 1/12. 

We carry out the CMB likelihood analysis for the mod- 
els g(Y) = c + c p YP with p = -0.5,0.5,1,2,3 by em- 
ploying the linear expansion given above. The observa- 
tional constraints shown in Fig. [5] with the bold lines 
are derived by using the theoretical values of n s , r, and 
rit given in Eqs. ([43]) and ([44]) with the relation ([79]) . 
In the (e, c 2 ) plane we also plot the border correspond- 



ing to the WMAP 7yr lower bound / 



cquil 
NL 



-214. The 



region above this border satisfies the observational con- 
straint of the non-Gaussianity. From Fig. [S]we find that 
there is no viable parameter space for p > 3 satisfying all 
the observational constraints. As long as the bold lines 
plotted in Fig. [S] are above the border corresponding to 
/S u = —214, the models with p < 3 can be compatible 
with the observational data. If the future observations 
can place the bound on /J^l" larger than —80, the mod- 
els with p > 1/2 can be ruled out (see Fig. [5]). 

Let us also discuss the case in which the function giY) 
is the sum of different powers of p. For example we con- 
sider the model 



g{Y) = co + Cl Y + c^Y- 1 



(80) 



This corresponds to the dilatonic ghost condensate in 
the presence of the exponential potential, i.e. P = 



9 



0.006 



0.005 



'» 0.003 




B. Class (ii) 

In the DBI model, inflation occurs in the ultra- 
relativistic regime (Y ~ 1/2) under the condition (|66[) . In 
this case the denominator of e in Eq. (|65[) is much larger 
than its numerator. Since the linear expansion around 
Y = 1/2 is not possible in such cases, we need to treat 
this class of models separately. Let us take the function 
of the form 



9<Y) 



Y 



[i + /on] , 



in which case Eqs. (|4"Tj) and (|4"2"j) give 

3Yf'(Y) 



1 + f(Y) - 2Yf'(Y) 
" s f'(Y) + 2Yf"(Y) ' 



(84) 

(85) 
(86) 



Figure 5: The thin lines show the relations between e and <? a 
for the models g(Y) = cq+c p Y p with (i) p — —0.5, (ii) p = 0.5, 
(iii) p = 1, (iv) p — 2, and (v) p = 3, in the regime e <C 1. 
The bold lines correspond to the observational constraints 
(95 % CL) on each model derived from the joint data analysis 
of WMAP 7yr, BAO, and HST. We also plot the boundary 
coming from > —214 as well as the curves corresponding 



to /N q L uil = -50 and /, 



cquil 
NL 



= -80. 



c Q X + cie^X 2 + C-ie"^. Substituting Eq. flSO]) into 
Eq. ([TO]), we obtain Y = -c /(2ci). The conditions (JTSJ) 
and ([76]) translate into ci(l — 4cic_i/cg) > and ci > 0, 
respectively. Since Yq = — cq/(2ci) > 0, we require that 



c <0, Ci>0, 4cic_i/cq < 1 



From Eq. (|T8|> we have 



c 



1 

12 



1 - 



4ci c_ 



(81) 



(82) 



If c_i > 0, then the tangent of the line (|82p gets smaller 
relative to that in the ghost condensate model. Figure 
[5] shows that the allowed parameter space tends to be 
narrower for smaller c 2 /e. The existence of a viable pa- 
rameter demands the following condition 



r 2 



< 0.1 



(83) 



The effect of the negative exponential potential V = 
— c_ie~ A ^ (with c_i > 0) needs to be suppressed to 



be consistent with the bound ([8"3]) . In contrast, the 
tangent of the line (|82|) gets larger than 1/12 when 
c_i < 0. The effect of the positive exponential poten- 
tial V = — c_ie ^ (with c_i < 0) makes it easier to 
satisfy the observational constraints. 



In the DBI model with f(Y) = y/1 — 2Y /c, we can realize 
inflation in the regime |/'(V)| = l/[|c|Vl - 2Y] < 1 and 
f(Y) < 1 with Y w 1/2, so that e < 1. 
We study the models ([54]) with 



From Eqs. 



/(r) = (i-2rr/ c . 

and ([55]) we have 
6mY 

c(l - 2Y) 1 -" 1 + 1 + 2(2m - 1)Y 
1 - 2F 



(87) 



1 - 2(2m- l)Y 



(89) 



We consider the case in which inflation occurs for the 
values of Y slightly smaller than 1/2, while satisfying 
the condition c(l — 2Y) 1 ~ m ^> 1. Since we require e ~ 
3m/[c(l-2r) 1 - m ] > and c 2 s ~ (l-2F)/[2(l-m)] > 0, 
we have either < m < 1 with c > or m < with c < 0. 
The following relation also holds between c s and e: 



„2(i-m) 



3m 



1 



c[2{l-m)Y- 



(90) 



ycquil 



In Fig. [BJwe plot the curve (|90l) in the (e, c^) plane for 
four different values of m (= 0.3,0.5,0.7,0.9). The left 
and right panels correspond to the cases c = 10 3 and c = 
10 4 , respectively. We also show the observational bounds 
constrained by n s , r, and nt (plotted as the bold lines) 
as well as the curves corresponding to — —214 and 

-50. 

When c = 10 3 there exists some allowed parameter 
space for the models with m < 0.7 (including the DBI 
model with m = 0.5), but for m > 0.9 the WMAP 7yr 
bound of the non-Gaussianity excludes the parameter re- 
gion constrained by the linear perturbations. For larger 
c the theoretical curves in Fig. [5] shift to the regions 
with smaller c s , so that the constraint from the non- 
Gaussianity tends to be more important. For c = 10 4 



10 



0.1 



0.01 



0.001 



0.0001 




0.01 



0.001 



a uil > "214 



0.0001 



1e-05 




Figure 6: The thin curves show the relation between e and c? a for the models (|84|l with (|87[l . The left and right panels correspond 
to c = 10 3 and c = 10 4 , respectively, with (i) m = 0.3, (ii) m = 0.5, (iii) m = 0.7, and (iv) m = 0.9. The bold curves represent 
the observational constraints (95 % CL) derived from the CMB likelihood analysis in terms of n s , r, and n t . We also plot the 



borders corresponding to the WMAP 7yr bound /] 



cquil 



-214 as well as f^ 1 



-50. 



the right panel of Fig. [5] shows that, the models with 
to > 0.5 do not have the viable parameter space satisfy- 
ing all the current observational constraints. If the future 
observations can reach the level of the lower limit of the 
non-Gaussianity with Uni!"'! = C(10), then it is possi- 
ble to place tighter constraints further (see the curves in 
Fig. [6] corresponding to = —50). 



VI. CONCLUSIONS 

We have studied the observational constraints on as- 
sisted k-inflation models in which the multiple scalar 
fields join an effective single-field attractor described by 
the Lagrangian P = Xg{Y) with Y = Xe X( ^. The canon- 
ical field with the exponential potential, P = X — ce~ x< ^ 
(i.e. g(Y) = 1 — c/Y), is one of the simplest ex- 
amples giving rise to assisted inflation. The effective 
slope A along the inflationary attractor is given by A = 

l/^f) i which is smaller than the slopes A; for 
each exponential potential. The same structure holds for 
the k- inflation models with the Lagrangian P = Xg(Y) 
for arbitrary functions of g{Y). 

Along the effective single-field attractor, the inflation- 
ary obscrvables are in general given by Eqs. (|43|) - (|45l) . 
In Sec. IIVI we have confronted four models of assisted 
inflation with the recent observations of CMB combined 
with BAO and HST. For the canonical field with the ex- 
ponential potential the effective slope A is constrained to 
be 0.086 < A < 0.228. The tachyon field needs to have 
a small kinetic energy relative to its potential energy for 
the realization of inflation, in which case the observa- 



tional bound on the variable Y is given by Eq. (|55|) . Since 
the field propagation speed c s is close to 1 in this case, 
the primordial non-Gaussianity remains to be small for 
the tachyon model. 

In the dilatonic ghost condensate model the non- 
Gaussianity provides additional constraints to those de- 
rived by the spectra of scalar and tensor perturbations. 
As we see in Fig. d the WMAP 7 yr limit /^ uiI > -214 
reduces the allowed parameter space of the parameter A. 
If the lower bound on f^ 11 reaches the level of —100 in 
future observations, it will be possible to rule out the dila- 
tonic ghost condensate model. In the DBI model the level 
of the non-Gaussianity depends on the field propagation 
speed c s as well as the constant c associated with the 
energy scale of the potential. For larger c, Ifrnf | tends 
to increase, so that the models can be constrained by the 
additional information coming from the non-Gaussianity. 
In fact the DBI model with c > 10 4 is excluded by the 
WMAP 7yr data. 

We have extended the analysis to more general func- 
tions g(Y) by classifying the assisted k- inflation models 
into two classes. The first class consists of the mod- 
els in which inflation occurs around Y = Yq satisfying 
the condition g(Yo) + Yog'(Yo) = 0. The representa- 
tive models of this class are g(Y) = cq + J2 P ^o c yY % ', 
which includes the dilatonic ghost condensate. From 
the CMB likelihood analysis combined with the non- 
Gaussianity bound we have found that the single-power 
models g(Y) = Co + c p Y p with p > 3 are ruled out. The 
second class consists of the models with the speed limit 
of the field, which includes the DBI model as a specific 
case. We have carried out the CMB likelihood analysis 
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for the functions g{Y) = -{c/Y)[l+{l-2Y) m / c] (m < 1) 
and showed that the models with larger m and c tend to 
be observationally disfavored by taking into account the 
non-Gaussianity bound. 

In this paper we have evaluated the inflationary ob- 
servables under the assumption that the solutions are 
on the assisted attractor described by the effective sin- 
gle field. In order to end inflation the solutions need to 
exit from this regime. This can be achieved by treating 
the validity of our Lagrangian P = Xg(Y) only within 
some limited range of field values. With some suitable 
modification of the Lagrangian it is possible to lead to 
the graceful exit of inflation [l3j . Another possibility is 
that k-inflation ends with a phase transition as in hy- 
brid inflation [2!| . It will be also of interest to study the 
case where the observed CMB anisotropics correspond 



to the epoch before the multiple fields join the inflation- 
ary attractor. In this case the trajectory in field space is 
curved, so that isocurvature perturbations can contribute 
to adiabatic perturbations [40[ . We leave these issues for 
future work. 
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